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Abstract 

We report I2 CO(7 =3-2) observations of 15 nearby elliptical galaxies carried out with the ASTE telescope. 
Thirteen were selected without regard to the presence of other tracers of cold interstellar matter. CO emission was 
detected from three of the galaxies, two of which are undetected by IRAS at 100 microns, suggesting that cold 
ISM may be present in more ellipticals than previously thought. The molecular gas masses range from 2.2 x 10 6 to 
4.3 x 10 8 Mq. The ratio of the CO(3-2) and (1-0) lines, R$i, has a lower value for elliptical galaxies than for spiral 
galaxies except for NGC 855, for which the value is close to the mean for spirals. The molecular gas in NGC 855 
has a mean density in the range 300 - 1000 cm~ 3 adopting a temperature range of 15 - 100 K. 
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1. Introduction 

Spiral galaxies contain substantial amounts of cold interstel- 
lar matter (ISM) and ongoing star formation, while cold ISM is 
but a minor component of elliptical galaxies. This is generally 
attributed to early exhaustion of the ISM in ellipticals during 
the early stages of galaxy formation. However, mass loss dur- 
ing stellar evolution (Faber & Gallagher 1976) can replenish 
the ISM to detectable levels in much less than a Hubble time, 
and ISM present in ellipticals can be attributed to stellar mass 
loss or to accretion from outside. 

The mass loss from evolved stars is estimated to be 1 .5M 
per year per 10 Lq of optical luminosity (Faber & Gallagher 
1976) and indeed is observed to take place as cold molecular 
gas (e.g., Teyssier et al. 2006). 

Thus, ISM contents comparable to those of spiral galax- 
ies could be accumulated in elliptical galaxies in of order 10 9 
years. However, Knapp et al. (1985) point out that the cold gas 
and stellar contents of elliptical galaxies are uncorrelated, and 
conclude that the cold gas comes from outside through events 
such as mergers. On the other hand, many elliptical galaxies 
contain of order 10 8 to 10 10 M P of hot gas, detected by its X- 
ray emission (Forman et al. 1985), which could originate in 
collisionally-heated mass shed by dying stars (e.g., Mathews 
& Brighenti 2003). 

Knowledge of the cold-gas content of normal elliptical 
galaxies is crucial to these questions. Several such surveys 
have been made (Lees et al. 1991; Sofue & Wakamatsu 1993; 
Wiklind et al. 1995; Knapp & Rupen 1996) using the CO(1-0) 
and CO(2-l) rotational lines as probes. However, all of these 
studies have concentrated on ellipticals which are detected at 
100 microns by IRAS, and are therefore biased. The present 
paper discusses observations of 15 nearby elliptical galaxies 
selected independently of their infrared emission. These ob- 
servations were made of the CO(3-2) line, which in addition to 
providing information on the presence of molecular gas allows 
an evaluation of its physical condition from multi-transition 



measurements (cf. Li et al. 1993; Vila-Vilaro et al. 2003). 
2. Observations 

We carried out 12 CO(/ =3-2) (hereafter CO(3-2)) ob- 
servations of fifteen elliptical galaxies using the Atacama 
Submillimeter Telescope Experiment (ASTE), a 10m antenna 
located at Pampa La Bola in the Atacama Desert in Chile at an 
altitude of 4800 m (Ezawa et al. 2004; Kohno 2005). The ob- 
servations were made on August 21-24 2006 from the remote 
ASTE operations room at the Nobeyama Radio Observatory 
(NRO) using the network operation system N-COSMOS 3 de- 
veloped at the National Astronomical Observatory of Japan 
(NAOJ) (Kamazaki et al. 2005). The front end is a cartridge- 
type 350 GHz receiver operated in double-sideband mode with 
an intermediate frequency (IF) between 5 and 7 GHz. The an- 
tenna temperature and sky extinction were measured using the 
standard chopper-wheel method. The antenna beam size is 22" 
at 345 GHz. We used three digital spectrometers with a band- 
width of 512 MHz (445 km s" 1 at 345 GHz) and 1024 chan- 
nels. The same center frequency was set for these spectrom- 
eters. A velocity resolution was 0.53 km s . The pointing 
was monitored by observations of the bright line sources AFGL 
3068 and o Ceti, and was accurate to 5". The main beam effi- 
ciency T|mb was measured each day by observations of W51D 
using the peak intensity measured by Wang et al. (1994) and 
was constant at 0.54. 

The observed galaxies were selected from the catalogue of 
Knapp et al. (1989) to lie in the right ascension range 22h - 7h 
and between declinations of —30 degrees — h20 degrees (so 
that the galaxies can be observed from both the Northern and 
Southern hemispheres), to have recession velocities less than 
5000 km s _1 and to be classified as ellipticals rather than as SOs 
or of later types. No other criterion was used for thirteen of the 
galaxies, so they are unbiased. In addition, two infrared-bright 
dwarf ellipticals, NGC 855 and NGC 2328, were observed. In 
these galaxies, CO(1-0) and CO(2-l) emission have been ob- 
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served and the line widths are less than 445 km s . The ob- 
served galaxies are listed in table 1 . 

The data reduction was carried out with tasks in NEWSTAR, 
a package developed at NRO. Individual bad data were flagged, 
the remaining spectra averaged, and first-order baselines were 
removed. The calibrated spectra were smoothed by 46 channels 
to a velocity resolution of 19.8 km s . The typical rms noise 
levels, ATmb, for the smoothed spectra were 8 mK. 

3. Results 

3.1. Spectra 

The spectra for all fifteen observed galaxies are shown in fig- 
ure 1 and the observational quantities are listed in table 2. We 
detected CO(3-2) emission from three galaxies, all for the first 
time (the observation of NGC 855 by Vila-Vilaro et al. (2003) 
is not sufficiently sensitive to detect emission at the level of that 
measured in this paper). 

The spectrum for NGC 1550 is highly asymmetric relative 
to the systemic velocity measured in optical observations. In 
order to examine whether this detection is real, we halved the 
original data into two independent groups, and then averaged 
the spectra belonging to each group. In the both cases, similar 
profiles were obtained as shown in figure 2. Hence, we can 
conclude that this detection is likely to be real. 

Although IRAS 100 micron flux in NGC 7464 is the high- 
est in our observed sample, CO emission has not been detected 
(Wiklind et al. 1995). NGC 7468 was detected in only CO(l- 

0) line at 22" beam (Wiklind et al. 1995), which is the same 
as our observing beam size. Since the CO(1-0) peak tempera- 
ture is about 20 mK, the CO(3-2) peak temperature would be 
about 4 mK, assuming the mean line ratio for elliptical galax- 
ies of /?3i = 0.21 (Vila-Vilaro et al. 2003). The rms noise of 7 
mK in our observation was not low enough to detect CO(3-2) 
emission. 

Although NGC 2328 was detected in both CO(1-0) and (2- 

1) lines (Wiklind et al. 1995; Lees et al. 1991), we could not 
detect CO(3-2) emission. The peak temperature of the CO(l- 
0) line was about 15 mK at 44" beam. If the size of the molec- 
ular distribution was less than 22", and if the Rt,\ was 0.21, the 
CO(3-2) peak temperature would be about 13 mK. The rms 
noise in our observation was too high (12 mK) to detect the 
line. 

Note that CO(3-2) emission could not be detected if a veloc- 
ity width was larger than the bandwidth of the spectrometers 
(445 km s- 1 ). 

3.2. Molecular Masses 

We assume the average ratio Rn of 0.21 to estimate the ex- 
pected CO(1-0) flux for two of the galaxies, and use the CO(l- 
0) data for NGC 855 from Wiklind et al. (1995). The molecular 
hydrogen mass is calculated from 

M^^ico^^eV. (i) 

4 K31 

The distance was calculated from the recession velocity V 
assuming a Hubble constant of 72 km s Mpc~' (Spergel 
et al. 2003). The conversion factor was taken to be Xqo — 
2.8 x 10 20 cm~ 2 K _1 km -1 s, an average of the values found 



for nearby Galactic molecular clouds by Arimoto et al. (1996). 
The resulting masses are listed in table 3 and cover a similar 
range to those measured for other elliptical galaxies (Lees et 
al. 1991; Sofue & Wakamatsu 1993; Wiklind et al. 1995). 

4. Discussion 

4.1. Line Ratios for NGC 855 

Both CO(1-0) and CO(2-l) emission have been detected 
from NGC 855 (Lees et al. 1991; Wiklind et al. 1995). 
Comparison with the observations of Wiklind et al. (1995) 
gives line ratios of Rt,\ = 0.5 and R21 = 0.8. The beam sizes 
used for the CO(3-2) and CO(1-0) are the same, 22", while the 
beam size for the CO(2-l) observation is 11". Therefore, the 
true value of R2\ might be a factor 1 to 4 lower, while Z?3i can 
be considered as the true value. 

The value of R31 is larger than the mean observed for 
three other elliptical galaxies (0.2, Vila-Vilaro et al. 2003), 
but similar to values observed for spiral and starburst galaxies 
(Mauersberger et al. 1999). 

4.1.1. Star Formation Efficiency 

The star formation efficiency (SFE) of NGC 855 is about 
8 x 10~ 10 yr" 1 (Vila-Vilaro et al. 2003), higher than that of 
most elliptical galaxies and similar to values observed for spiral 
galaxies (Komugiet al. 2005). This suggests that/?3i correlates 
with star formation efficiency, as noted from mapping of M83 
by Muraoka & Kohno (2006a) and Muraoka et al. (2006b). 

4.1.2. Modeling the Molecular ISM 

Since we have data for three CO rotational in NGC 855, the 
mean properties of the molecular ISM can be modelled. First, 
we assume LTE and that all three CO lines are optically thick. 
The brightness temperature, 7b, is then 

*-£(K£)-0~4-*(&H> 

where h is Planck's constant, v is frequency, k is Boltzmann's 
constant and T is the kinetic temperature. Applying equation 
(2) to the observed line ratios for NGC 855 implies a kinetic 
temperature of 8 K, which is far lower than the dust tempera- 
ture of 33 K estimated by Wiklind et al. (1995). 

Accordingly, we modelled the emission using the Large 
Velocity Gradient (LVG) model (Goldreich & Kwan 1974). We 
calculated line ratios for kinetic temperatures from 10 - 100 K, 
H2 densities from 10 to 6.3 x 10 6 cm~ 3 , and CO abundances 
XcoidV/dr)- 1 of 10~ 6 , 10~ 5 and 10~ 4 pc km" 1 s. The re- 
sulting line ratios Ri,\ and R21 are shown in figure 3. The data 
and models match well for T > 15 K, compatible with the dust 
temperature. The corresponding density is in the range of 300 
- 1000 cirT 3 . The CO(1-0) line is optically thick at this den- 
sity, as was assumed in estimating the molecular mass from 
Equation (1). 

4.2. Volume and Size of Molecular Distribution 

Using the calculated mean density, we can estimate the vol- 
ume of the molecular gas and its diameter, assuming a spherical 
distribution. We used a mean H2 density of 650 cm~ 3 for all 
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three detected elliptical galaxies. The results are given in table 
3. The molecular content in NGC 855 is as much as a single 
Galactic Giant Molecular Cloud (GMC), which can support the 
existence of the star formation in this galaxy (e.g., Vila-Vilaro 
et al. 2003). 

However, the size is smaller than the beam size of CO(2-l) 
(11"), which corresponds to 440 pc at 8.3 Mpc. Therefore, R21 
would be 4 times smaller (i.e., R 2 \ = 0.2), if NGC 855 had only 
one GMC. The LVG calculation cannot satisfy R$\ = 0.5 and 
^21 = 0.2 at the same time, if this was the case. The curve of 
R31 = 0.5 always coincides with that of ^21 ~ 0.8. 

Considering these matters, the molecular distribution is 
likely to be more extended than the beam size 22", forming 
smaller clumps rather than one GMC. 

4.3. Comparison with IRAS fluxes 

Two of the elliptical galaxies (NGC 990 and NGC 1550) in 
which we detected molecular gas are not detected by IRAS, 
though IRAS lOOjUin upper limit of NGC 990 is pretty high. 
Thus cold molecular ISM may be more common in ellipticals 
than previously thought. 
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5. Summary 

We achieved 12 CO(/ =3-2) line survey of 15 elliptical 
galaxies and detected confident CO emissions from three el- 
lipticals. The molecular gas masses were estimated to be 
2.2 x 10 6 -4.3 x 10 8 M o . 

Line ratio, Z?3i, of NGC 855 was calculated to be 0.5. 

We made line ratio analysis using R31 and R21 of NGC 855. 
The LVG model showed that H 2 density was 3 x 10 2 - 1 x 10 3 
cm~ 3 in the temperature range of 15 - 100 K. 

CO(7 =3-2) emission was detected even in elliptical with 
faint IRAS 100 /jm flux. More ellipticals might be detected in 
the CO lines than previously thought. 

We are grateful to all the ASTE staff for their dedicated sup- 
port for our observation. We thank the referee Dr. Knapp for 
substantial corrections of English as well as invaluable scien- 
tific comments, which greatly improved our paper. We thank 
Dr. T. Sawada and Dr. S. Takano for fruitful discussions 
and advices for the observations. This study was financially 
supported by MEXT Grant-in-Aid for Scientific Research on 
Priority Areas No. 15071202. 
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Table 1. Target list 



Name 


R.A.(B1950) 


Dec.(B1950) 


Morph. 


B7 1 


Vhel-opt 




h 


m 


s 


d 


m 


s 




mag. 


km s _1 


First samples 


NGC596 


01 


30 


22 


-07 


17 


18 


El 


11.8 


1817 


NGC636 


01 


36 


36 


-07 


45 


54 


E3 


12.4 


1805 


NGC821 


02 


05 


41 


+ 10 


45 


32 


E6 


11.7 


1716 


NGC990 


02 


33 


36 


+ 11 


25 


28 


E 


13.9 


3508 


NGC1453 


03 


43 


57 


-04 


07 


36 


E2 


12.6 


3906 


NGC1550 


04 


17 


02 


+02 


17 


25 


E 


14.0 


3689 


NGC1600 


04 


29 


12 


-05 


11 


30 


E4 


12.0 


4687 


NGC7458 


22 


58 


55 


+01 


29 


05 


E 


13.9 


4981 


NGC7464 


22 


59 


25 


+ 15 


42 


17 


Elp 




1877 


NGC7468 


23 


00 


30 


+ 16 


20 


08 


E3p 


14.0 


2089 


NGC7619 


23 


17 


43 


+07 


55 


57 


E2 


12.3 


3747 


NGC7626 


23 


18 


10 


+07 


56 


35 


E 


12.3 


3450 


NGC7785 


23 


52 


45 


+05 


38 


11 


E5 


12.6 


3826 


Second samples 


NGC855 


02 


11 


10 


+27 


38 


36 


E 


13.3 


600 


NGC2328 


07 


01 


01 


-41 


59 


42 


E/SO 


13.3 


1159 



These parameters are taken from Knapp et al. (1989). 
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Fig. 1, CO Spectra of all the observed elliptical galaxies. 
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Fig. 2. Spectra for NGC 1550. We halved the original data into two independent groups. Presented spectra are obtained by averaging the spectra belonging 
to the individual groups. 



Table 2. Observational quantities 



Name 


fl00/im 


Vco 


AV C0 


^peak 


Ico 


A7mb 


comment 




Jy 


km s 


km s _1 


mK 


Kkms- 1 


mK 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


First samples 


NGC596 


< 0.300 






8 




8 


no detection 


NGC636 


< 0.399 






12 




5 


no detection 


NGC821 


0.440 






19 




6 


no detection 


NGC990 


<2.24 


3531 


40 


16 


0.59 


3 


detection 


NGC 1453 


0.670 






12 




7 


no detection 


NGC 1550 


< 0.651 


3631 


198 


37 


4.9 


5 


detection 


NGC 1600 


0.170 






26 




15 


no detection 


NGC7458 


< 0.192 






12 




8 


no detection 


NGC7464 


6.71 






16 




8 


no detection 


NGC7468 


1.65 






7 




7 


no detection 


NGC7619 


0.630 






13 




11 


no detection 


NGC7626 


< 0.339 






17 




7 


no detection 


NGC7785 


< 0.264 






20 




11 


no detection 


Second samples 


NGC855 


3.27 


613 


20 


17 


0.36 


4 


detection 


NGC2328 


3.77 






25 




12 


no detection 



Col. 1 . — Galaxy name. 

Col. 2. — IRAS 100/jm flux. The upper limit was taken to be 3a. 
Col. 3. — Central velocity. 

Col. 4. — Velocity width (Full width at zero intensity; FWZI). 

Col. 5. — Peak brightness temperature. 

Col. 6. — Integrated intensity. 

Col. 7. — Rms noise. 

Col. 8. — Comment on detection. 



Table 3. Inferred Properties of Molecular Gas 



Name 


M m 


Volume 


Diameter 




M 


3 

pc 


pc 


NGC990 


2.7 x 10 8 


(8.3 x 10 b ) 


(250) 


NGC 1550 


4.3 x 10 8 


(1.3 x 10 7 ) 


(300) 


NGC855 


2.2 x 10 6 


6.9 x 10 4 


(51) 
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Fig. 3. Results of LVG calculations adopting CO abundance of Xco/ (dV/dr) = 10 5 pc km 1 s. The solid and dotted lines denote R31 and R21, 
respectively. The thick solid and dotted lines denote R3i=0.5 and i?2i=0.8, respectively. 



